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Abstract. Since the discovery of brown dwarfs in 1994, and the discovery of dust 
cloud formation in the latest Very Low Mass Stars (VLMs) and Brown Dwarfs (BDs) 
in 1996, the most important challenge in modeling their atmospheres as become the 
understanding of cloud formation and advective mixing. For this purpose, we have de- 
veloped radiation hydrodynamic 2D model atmosphere simulations to study the forma- 
tion of forsterite dust in presence of advection, condensation, and sedimentation across 
the M-L-T VLMs to BDs sequence (T eS = 2800 K to 900 K, Freytag et al. 2010). We 
discovered the formation of gravity waves as a driving mechanism for the formation 
of clouds in these atmospheres, and derived a rule for the velocity field versus atmo- 
spheric depth and T e s, which is relatively insensitive to gravity. This rule has been used 
in the construction of the new model atmosphere grid, BT-Settl, to determine the micro- 
turbulence velocity, the diffusion coefficient, and the advective mixing of molecules as 
a function of depth. This new model grid of atmospheres and synthetic spectra has 
been computed for 100,000K > T eff > 400 K, 5.5 > logg > -0.5 and [M/H]= +0.5 to 
-1.5, and the reference solar abundances of Asplun det al.l ([2009). We found that the 
new solar abundances allow an improved (close to perfect) reproduction of the photo- 
metric and spectroscopic VLMs properties, and, for the first time, a smooth transition 
between stellar and substellar regimes — unlike the transition between the NextGen 
models from Hauschildt et al. 1999a,b, and the AMES-Dusty models from Allard et 
al. 2001). In the BDs regime, the BT-Settl models propose an improved explanation 
for the M-L-T spectral transition. In this paper, we therefore present the new BT-Settl 
model atmosphere grid, which explains the entire transition from the stellar to planetary 
mass regimes. 



The Impact of the new Solar Abundances on VLMs spectral properties 



The modeling of the atmospheres of very low mass stars (hereafter VLMs) has evolved 
with the development of computing capacities fro m an analytica l treatment of the trans- 
fer equation using moments of the r adiation field (|Allardlll990l), to a line-by-lin e opac- 
ity sampling in spherical sym metry (lAUard etalJI 19971 : lHauschildt et alJ[l999 a.b). and 
finally 3D radiation transfer dSeelmann et all [2010 ). In parallel to detailed radiative 
transfer in an assumed static environment, hydrodynamical simulations have been de- 



1 



2 



Allard, E, Freytag, B., and Homeier, D. 



veloped to reach a realistic representation of the granulation and the line profiles shifted 
and shaped by the hydrodynamical flow of the sun and sun-like stars (see for details the 
review in a special isssue of the Journal of Computational Physics by Freytag et al. 
2011) by using a non-grey (multi-group binning of opacities) radiative transfer using a 
pure blackbody source function (scattering is neglected). 



o 




Ludwig 



_i i i_ 



_i i i_ 



1.5 



2.5 



Wavelength (ytxm) 



Figure 1 . Model atmosphere synthetic spectrum ver sus VB10 IR SEP, using dif- 
ferent water opacity profiles over the years (lLudwidll97l lj^ gensen et al . 2001; 
ISchrvber et al.lll995l:IPartridge & Schwenkdl 19971) . 



The model atmospheres and synthetic spectra have also been made possible thanks 
to the development of realistic opacities calculated often ab initio for the needs of an 
accurate account of their cooling and heating effects in the internal atmospheric lay- 
ers where temperatures close to 3000 K can prevail. For some time, the remaining 
discrepancies in the model synthetic spectra were believed to be due to incomplete wa- 
ter vapor line lists, both in temperature and in the rotation quantum number / of the 
molecular simulations. In FigUl model atmosphere synthetic spectra are shown using 
differe nt published w ater vapor opacity profil es such as hot flames laboratory experi- 
ments (lLudwiglll97lh . empirical calculations ( J0rg ensen et al.ll200lh . and ab initio cal- 
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culations (ISchryber et al.lll995l : IPartridge & Schwenkdfl997h based on independently 
measured interaction potential surface. As can be seen from FigH] where the models 
are compared to the infrared spectrum of an M8 dwarf (VB10), the water vapor opac- 
ity profile which shape this part of the spectrum has strongly change over time with 
the improvement of computational capacities and a better knowledge of the interaction 
potential surface. But in general, all opacity profiles converge in predicting an over- 
opacity (or lack of flux in the model) in the K bandpass. The UCL opacity profile — 
likely because of its incompleteness — could allow a seemingly correct J - K color 
while not allowing a detailed spectral comparison of the NextGen models (Allar detall 
119941). It became however clear that the more rec ent versions of the water vapor profile 
( Partridge & Schwenkdfl997l : iBarber et al.ll2006h all agree in establishing a systematic 
lack o f flux in the K bandpass in the models. This is t he reason why thelAUard et aU 



(1200 lh AMES-Cond/Dusty model grids, based on the IPartridge & Schwenkd (1199^/ 
water vapor opacity profile were never proposed for the study of VLMs, and reserved 
for the study of the limiting properties of brown dwarfs. 
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Figure 2. A BT-Settl synthetic sp ectrum with TefF=2900K, logg=5.0, and solar 
metallicity by Asplun d"et all d2QQ9h is shown in yellow compared to the SED of 
the red dwarf GJ 866 (cyan curve). Both are plotted in absolute flux, i.e no flux 
adjustment at a specific wavelength is needed. The infrared SED (top panel) is now 
perfectly reproduced by the model. The agreement is particularly good in the Wing 
Ford FeH bands near 0.99 yum while some discrepancy prevail in the missing CaOH 
bands around 6000A. Telluric absorption is corrected from the observations in red. 



In this regard, the new solar a bundances based on radiation hydrod ynamical sim- 
ulations of the solar atmosphere (lAsplund et al.ll2009l : ICaffau et al.ll2010 ) help as they 
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predict an oxygen a bundance of 0.3 dex (a factor of 2) lower than the previously used 
solar abundances of iGrevesse et alJ (Il993h . The results are shown in FigEl where the 
new BT-Settl model atmos phere synthetic spectr a for a T q q = 2900 K, logg=5.0 and 
solar abundances according lAsplund et alJ (l2009h are compared to the spectrum of an 
M6 dwarfs, GJ866 (kindly provided by M. Bessel, Mt-Stromlo Obs.). For the first time, 
we find a perfectly fitting spectra distribution across the near-IR to infrared spectral re- 
gion (the model is the yellow line). The agreement is also excellent in the optical to red 
part of the spectrum in particularly in the FeH Wing Ford bands near 0.99 yum, and in 
the VO bands thanks to line lists provided by B. Plez (GRAAL, Montpellier, France). 
Missin g opacities are, how ever, affecting still the spectral distribution (CaOH bands), 
and the Allard et al. (2000) TiO line list becomes now too corse compared to progress 
made with water vapor. 
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Figure 3. Estimated Teif f or M dwarfs by Casagra nde et all (l2008[) and brown 
dwarfs by iGolimowski et all (|2004|) are reported as a function of J-K short. Q ver- 
plotted are the NextGen model isochrones for 5 Gyrs (iBaraffe et al .11 19971 1 19981) us- 
ing varius generations of model atmospheres, starting with the Ne xtGen (black line) , 
pursuing with the limiting case AMES-Cond/Dusty grids by Alla rd et all ([2001) 
(blue and red line respectively), and finishing with the BT-Settl models using the 
lAsplund et al.l (|2009|) solar abundances (green line). 



One can see from Fig J3] that the NextGen models systematically and increasingly 
overestimate T G s through the lower main sequence, while the AMES-Cond/Dusty 
models were underestim ating 7W compared to t he averaged empirical determinations 
of r e ff of i ndividual stars ( Casagrande et alJ[2008h . This situation is relieved when using 
the newer LAsplund et aD (|2009) abundances, and the BT-Settl models now agree fairly 
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well with most of the empirical estimations of T q q. Evolution models are currently 
being prepared using the BT-Settl model atmosphere grid. 



2. Dust formation in late type VLMs and Brown Dwarf Atmospheres 

One of the most important challenge in modeling the atmospheres and spectral proper- 
ties of VLMs and brown dwarfs is the formation of dust clouds and its associated green- 
house effects making the infrared colors of late M and early L dwarfs extremely red 
compared to colors of low mass stars. The cloud composition, according to equilibrium 
chemistry, is going from zirconium oxide (ZrC>2), to refractory ceramics (perovskite 
and corundum; CaTi03, AI2O3), to silicates (forsterite; Mg2SiC>4), to salts (CsCl, RbCl, 
NaCl), and finally to ices (H2O, NH3, NH4SH) a s brown dwarfs cools down with time 
from M through L, T spectral types and beyond (Allar( ietaDl200ll : lLodders & Fegley 



2006). Many cloud models ha ve been constructed to address this problem in brown 
dwarfs over the past decade (see lHelling et al.l l2008. for a review on the subject). How- 



ever, none treated the mixing properties of the atmosphere, and the resulting diffusion 
mechanism realistically enough to reproduce the properties of the spectral transition 
fro m M through L and T spect ral types without changing cloud parameters (for exam- 
ple lAckerman&MM^I |2QQl|). It is in this context that we have decided to address 



the issue of mixing and diffusion by 2D Radiation HydroDymanic (hereafter RHD) 
simulations of VLMs and brown dwarfs atmospheres, using the Phoenix opacitie s in a 
multi-group binning, and forsterite geometric cross-sections dFreytag et al.1 120101) . We 
found that gravity waves have a decisive role in clouds formation in brown dwarfs, 
while around T e ff < 2200 K the cloud layers become optically thick enough to initiate 
cloud convection, which participate in the global mixing. Overshoot is also important 
in the mixing of the largest dust particles (see paper by D. Homeier in this book). In 
Figgl 3D RHD simulations (a 350 x 350 x 170 km 3 box at the surface of the star) are 
shown for dwarfs with T q q = 2600 K, 2200 K, and 1500 K from top to bottom. For each 
simulation, two snapshots are shown side-by-side to illustrate the intensity variation 
due to cloud formation and granulation. The 2600 K case shows no or negligible dust 
formation, while dust formation progresses to reach optically thick density at around 
2200K, before sedimenting out again towards the 1300K regime. The T e ff =1500 K 
case illustrate the importance of gravity waves, where the minima of the waves reach 
condensation levels while the maxima remain in condensed phase. The box of simula- 
tions are too small compared to the radius of the star to show adequately the variability 
of these objects, however. See iFreytag et al.1 ([2010) and B. Freytag's poster paper in 



these proceedings for details. 

These simulations permitted to account in Phoenix for the advective forces bring- 
ing fresh condensible material from the hotter lower layers to the cloud forming layers. 
Our cloud model is based on the condensation and sed imentation tim escales from a 
study of Earth, Venus, Mars and Jupiter atmospheres by iRossowl (Il978h . However we 
had to compute the supersaturation pressure from our pre-tabulated equilibrium chem- 
istry in order to obtain the correct amount of dust formation (as opposed to use an 
approximate value cited in Rossow, 1978). We then solved the cloud model and equi- 
librium chemistry in turn layer by layer inside out to account for the sequence of grain 
species formation as a function of cooling of the gas. One can see from Figj5]that the 
late-type M and early type L dwarfs behave as if dust is formed nearly in equilibrium 
with the gas phase with extremely red colors in some agreement with the BT-Dusty 
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Figure 4. 3D HRD simulations using C05BOLD dFrevtag et al.ll201Qh of a 350 x 
350 x 170 km 3 of atmosphere at the surface of, from top to bottom, 2600K, 2200K 
and 1500K VLMs and brown dwarfs of logg=5.0, and solar metallicity. For each 
model two intensity snapshots are shown in time to illustrate the intensity variability. 



models. At the low TefF regime dominated by T-type dwarfs, the AMES-Cond mod- 
els also appear to provide a good limitation for the brown dwarf colors. However, the 
chemistry of T dwarfs can be very far from the equilibrium in the models. The BT-Settl 
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Figure 5 . Same plot as Figj3]but zooming out and extending into the brown dwarfs 
region of the diagram. This region below 2500K is dominated by dust formation 
(essentially forsterite and other silicates). The limiting case AMES-Cond/Dusty 
models atmosphere provide a description of the span in colors of the brown dwarfs 
in this diagram. 



models, which account for a cloud model, dynamical mixing from RHD simulations, 
a s upersaturation comput ed from pre-tabulated equilibrium chemistry calculations, and 
the lAsplund et alJ (l2009h solar abundances, manage to reproduce the main sequence 
down to the L-type brown dwarf regime, before turning to the blue in the late-L and 
T dwarf regime. The models used an age of 5 Gyrs, and in the case of the BT-Settl 
models, a younger age of a few Gyrs would easily reproduce the reddest brown dwarfs. 



3. Summary and Futur Prospects 



We propose, in this paper, a new model atmosphere grid, named BT-Settl, c omputed us- 



ing th e atmosphere code Phoenix which has bee n updated, compared to the lAllard et al. 



(12001 ) AMES-Cond/Dusty models, for: i) the Barber et alJ (12 006) BT2 water opac- 
ity line list, ii) the solar abundances revised by Asplun d et al. (2009), and iii) a cloud 
model accounting for supersaturation and RHD mixing. It is covering the whole range 
of VLMs and brown dwarfs and beyond: 1000,000 K < r eff < 400 K; -0.5 < logg < 
5.5; and +0.5 < [M/H] < -4.0, including various values of the alpha element enhance- 
ment. Only the confrontation of the models using spectral synthesis will allow to define 
the content of oxygen and alpha elements of VLMs and brown dwarfs. But it is clear 
that these objects are excellent constraint for the solar oxygen abundance. The models 
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are available at the Phoenix simulator website "http://phoenix.ens-lyon.fr/simulator/" 
and are in preparation for publication. However, the interior and evolution models are 
expected for the second half of 2011. 

In order to say something about the spectral variability of VLMs and brown dwarfs, 
3D RHD simulations of "the star in the box" with rotation will be required. This is our 
current project supported by the French "Agence Nationale de la Recherche" for the pe- 
riod 2010-201 5. Rotation is already m odeled for a scaled down model of the Sun using 
C05BOLD (ISteffen & Frevta g 2007) and can be applied to brown dwarf simulations. 
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